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'nI" . Abstract 

o 

■"nI" ' We prove that the Gross-Pitaevskii equation correctly describes the ground state 

i energy and corresponding one-particle density matrix of rotating, dilute, trapped Bose 

I gases with repulsive two-body interactions. We also show that there is 100% Bose- 

' Einstein condensation. While a proof that the GP equation correctly describes non- 

^ |. rotating or slowly rotating gases was known for some time, the rapidly rotating case was 

unclear because the Bose (i.e., symmetric) ground state is not the lowest eigenstate of 
the Hamiltonian in this case. We have been able to overcome this difficulty with the 
aid of coherent states. Our proof also conceptually simplifies the previous proof for the 
slowly rotating case. In the case of axially symmetric traps, our results show that the 
appearance of quantized vortices causes spontaneous symmetry breaking in the ground 
^ , state. 

Dedicated to Jakob Yngvason on the occasion of his 60^'^ birthday 



1 Introduction 

In this paper we sliow that a dilute, rotating Bose gas is correctly described by the Gross- 
Pitaevskii (GP) equation. We also show that there is 100% Bose-Einstein condensation 
(BEC) into a solution of the GP equation in a suitable limit. These conclusions were hereto- 
fore unproved, and it might not be an exaggeration to say that they were even conjectural, 
primarily because of the unusual situation (proved in [25] ) that the absolute ground state of 
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the Schrodinger Hamiltonian is not the bosonic ground state in the rapidly rotating case, as it 
is in the case when there is Httle or no rotation. In other words, the vortices seen in rotating 
gases are not properties of the absolute ground state but are, instead, true manifestations 
of the bosonic symmetry requirement. 

If the GP equation correctly describes the physics of a rotating gas (as we show here), 
then it also shows the superfluidity of such a gas, as will be discussed below. In the case of a 
cylindrically symmetric trap potential, the rotational symmetry is broken when more than 
one vortex is present; the GP equation must describe this broken rotational symmetry and, 
therefore, it must have multiple minimum energy solutions in this case. 

The key mathematical tool employed here is coherent states. Our work is based on the 
results of JS| and the observation there that one can make a c-number substitution for many 
boson modes (not just one, as in Bogoliubov's method) without significant error provided 
the number of such modes is of lower order than N, the number of particles. 

As in our previous work ^1 El UHl on dilute, trapped Bose gases we start with 
the Hamiltonian for bosons 

N 

1=1 l<i<j<N 

where Hq is the one-body part of the Hamiltonian and is the two-body repulsive interac- 
tion. These terms and the GP limit are described as follows. 

1. The GP limit: We want to fix the external trapping potential but let tend to 
infinity. To retain the notion of a dilute gas in this situation we let the interparticle potential 
depend on A^ in such a way that a^^ , the two-body scattering length of , is related to A^ 
by the condition that 

Noj^ = a is fixed. (2) 

In this limit the three components of the energy (kinetic, trapping potential and interaction 
potential) scale in the same way and are all of the same order of magnitude. We call this 
the GP limit. It is this limit that will lead to the GP equation 0. 

2. The two-body potential: We choose a radial two-body potential w{x) such that 
w{x) > (this is an important restriction for our methods) and such that w{x) = for 
l^l > Rq (this finite range condition is a technical restriction for simplicity and can be relaxed 
if need be). We note that integrability of w{x) is not assumed here, w{x) is even allowed to 
have a hard core. The scattering length of w is a (i.e., the solution to [—^A-'rw{x)]f{x) = 
with /(oo) = 1 satisfies f = 1 — a/\x\ for |x| > Rq). The actual two-body potential in (P), 
given by 

v^{x) = N'^w{Nx) , (3) 

has scattering length = a/N . 

3. The one-body Hamiltonian: We work, as usual, in the rotating coordinate system, 
in which case the kinetic energy has to be supplemented by a term —Q. ■ (pAx) = p - Ax), 
where is the angular velocity vector, and p = —ikV. It is convenient to add and subtract 



EHLRS November 22, 2005- GP-rot 



3 



a term ^(ri A x)^ and thereby write 

Ho = ^{p + A{x)f + V{x) (4) 

with A{x) = rriQAx. Then V is the trapping potential (which might or might not have some 
geometric symmetry) minus A x)"^. It is well known that we must have V{x) — > oo as 
|x| oo, for otherwise the system will fly apart. We can also assume that V > without 
loss of generality. Actually, for technical reasons we require just a little more, namely 
V{x) > Ci ln(|x|) — C2 for some positive constants Ci and C2. (This condition can probably 
be relaxed a bit. What we actually need is that Tr e"(^~^(^)) and Tr |^(x)|^e"(^~^('')) are 
finite for a large enough, for some s > 2. We will show in the appendix that this is fulfilled 
under the stated assumption on V.) 

We note that in the rotating coordinate system, the velocity at x is not p/m but rather 
V = ih~^[HQ, x] = p/m+VLf\x. The angular velocity around the axis is r2-(t)Ax)|x-'-|~^|f]|~^, 
where \x'^\ is the distance to the axis. In a cylindrically symmetric state if) we have 
• (p A x) ^ = and, therefore, the angular velocity is 17, not zero. In the fixed frame the 
angular velocity is 17 — = 0. In other words, the system in such a state is not rotating. As 
long as ri is small enough, the GP ground state is cylindrically symmetric and hence there 
is no rotation; this is a manifestation of superfluidity. In order to have rotation at least one 
vortex must form. This is a typical property of superfluids. 

Henceforth, we use units in which h = 2m = 1. We also note that the modiflcation of 
the kinetic energy in @ is mathematically just like that caused by a uniform magnetic fleld 
with vector potential A (and e/c = 1). There is nothing special about A{x) = mil A x as 
far as the mathematics is concerned, so one could have an arbitrary A without disturbing 
our analysis, provided it did not grow too fast at inflnity. One could think, for example, of 
applying a magnetic fleld to the system, but then our particles would have to be charged 
and the attendant Coulomb interaction would nullify the treatment of the system as a dilute 
gas with short range interaction. On the other hand, we could allow our particles to have 
a magnetic moment ("bosons with spin") and our analysis would easily extend to this case. 
The ground state energy depends in a non-trivial way on the total spin when there is rotation 
j22j . even in the absence of a magnetic fleld. This is due to the symmetry requirement of 
the wave function, whereby the symmetry of the spin part determines the spatial symmetry 
(see, e.g., [S]). We will not pursue this topic further in this paper. 

Our analysis is carried out here for three-dimensional gas particles, but the same ideas 
apply to a two-dimensional gas. There will be changes, of course, because the notion of 
scattering length is different in 2D and because the energy per particle of a homogeneous 
gas of low density g is not AttqUs as in 3D but rather 47r£)/| log Qa1\. (Here, is the unsealed 
scattering length of the interaction potential, which is held flxed in the thermodynamic limit 
for the homogeneous gas.) Thus, the GP equation will be a little different, but the conclusion 
will be the same: The only effect of rotation is to replace p'^ by \p + in the GP equation 
derived in ^S]. In order to keep this paper manageable we do not discuss the 2D case, but 
the interested reader can easily combine the results in |23j and the present paper. 
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The Hamiltonian Hj^ acts on L^(M^^) but we are interested in its restriction to the 
bosonic subspace of L^(M^^), namely to permutation symmetric functions. We denote the 
ground state energy of in the bosonic sector by Eq{N), and we keep in mind that this 
might be larger than the absolute ground state energy of when no permutation symmetry 
is imposed. 

We turn now to the GP equation, which originates from the GP energy functional for a 
complex- valued function of one variable x G R^. For a > 0, the GP energy functional is 
given by 

f = {mo\<P) + 4vra / |0(x)|^dx . (5) 



It can easily be shown ^3] that S*^^ [0] has a minimum over all (j) with \\4>\\2 = 1 and this mini- 
mum energy is denoted by E'^^{a). (We use the standard notation \\(j)\\p = [f \4>{x)\p dxj .) 
There might be several minimizers (and there surely will be when the trap has axial sym- 
metry and a is large |221 123j ) but each minimizing (p will satisfy the GP equation 

{-iV + A{x)f(l){x) + V{x)(t){x) + 87ra\(l){x)f(t){x) = ficpix) , (6) 

where is the chemical potential (i.e., the energy per particle to add a small number of 
particles). Note that = £''^^(0) + Ana J \(j){x)\^dx > E'^^{a) because of the quartic non- 
linearity. 

Our main theorem concerning the bosonic ground state energy of is the following. 
Theorem 1. With a denoting the scattering length of w, we have 

hm Mn = ^CP(,). 

A"— >oo iV 

In 1231 it was shown that 



limsup^<i?GP(a), (8) 

and, therefore, it remains only to prove a lower bound to liminfjv^oo Eq(N)/N of the right 
form, which we do here. 

The GP energy minimizer(s) (p also tells us something about the density (diagonal and 
off-diagonal) and about Bose-Einstein condensation in the ground state of H^^, or any 
approximate ground state. We call a sequence of bosonic A^-particle density matrices 7^ 
an approximate ground state if limjv^oo ^~^Tr = E'^^{a). The reduced one-particle 

density matrix of 7^^ will be denoted by 7^^) . 

We would like to suppose that, as — > 00, 7^-^^ converges to some 7 and that 7 = \4>){4>\, 
where (/> is a solution to the GP equation. This would be 100% Bose-Einstein condensation 
into the GP state and was proved to occur in the non-rotating case ,1%. The difficulty in 
the rotating case is that the solution to the GP equation might not be unique (as it is in the 
non-rotating case), in which case the limit 7 need not be a pure state. We would expect, 
however, that 7 is always a convex combination of pure GP states, i.e., 7 = X^j '^i I 
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where (pi is a solution to the GP equation and Aj = 1. (This, of course, is not the same 
as the much weaker and less interesting statement that 7 is a convex combination of terms 
of the form 1-0) (-01, in which -0 is a linear combination of GP solutions instead of being equal 
to just one GP solution.) Unfortunately, as in the case of a cylindrically symmetric trap, 
the set of GP states might not be countable, and so the summation J2i must be replaced 
by some kind of integral. This accounts for the rather abstract Theorem |2l below. In any 
event, this theorem tells us that there is always 100% condensation, even if the system has 
a wide choice of states into which to condense. 

Note that 7^-^) is a positive trace class operator on the one-particle space L^(]R^), and 
we choose the normalization Tr7^^) = 1 for convenience. (The conventional normalization 
is Tr7^^'* = A^.) By the Banach-Alaoglu Theorem, any sequence 7^^-* will have a subse- 
quence that converges to some 7 in the weak-* topology, i.e., limAr_>oo Tr ^7^^^ = Trj47 
for all compact operators A. This convergence will even hold in the norm topology, i.e., 
limAT^oo Tr [7^^) — 7I = by compactness. More precisely, since the 7^^) are the one-particle 
density matrices of approximate ground states, we have (using the positivity of the interac- 
tion potential in Hj^) Tr//o7j^^^ < const, independently of A^. Hence also \/ Ho'jj^^ ^'TIq 
\/ Ho'jy/ Hq in weak-* sense, i.e., Tr A-^ Ho'yj^^ V Hq — > Ti A^/H()^^/HQ for all compact A. 
Since Hq^ is a compact operator, this implies that Tr7^^-' ^ Tr7 as — > 00 (simply use 
A = Hq^ above). For positive operators, weak-* convergence plus convergence of the trace 
implies norm-convergence |271 125j . 

We denote by T the set of all 7's that are limit points of one-particle density matrices of 
approximate minimizers. That is, 

r = <7 : there is a sequence 7^^, lim —TrHj^jj^ = E'^^{a), lim 7^^-' = 7I . (9) 

I N^oo N N^oo J 

As remarked above, the convergence 7^^^ — > 7 can either mean weak-* convergence or norm 
convergence. Note that, in particular, norm convergence implies that Tr7 = 1 for all 7 G P. 

Theorem 2. The set T of one-particle density matrices of approximate ground states, as 
defined in 0), has the following properties. 

(i) T is a compact and convex subset of the set of all trace class operators. 

(a) Let Text C P denote the set of extreme points in T. (An element 7 G P is extreme if ^ 
cannot be written as 'j = 07^^ + (1 — a)^^ with 7^ 2 G P, 7^ 7^ 72; ^^^^ < a < 1.) We 
have Pext = {|0)(0| '■ S'^^i'p] = E'^^{a)}, i.e., the extreme points in P are given by the 
rank-one projections onto GP minimizers. 

(Hi) For each 7 G P, there is a positive (regular Borel) measure d/i^, supported in Pext, with 
Jp ^ dfi-y^cp) = 1, such that 

7 = / d/.^(0)|</.)(0|, (10) 

Text 

where the integral is understood in the weak sense. That is, every 7 G P is a convex 
combination of rank-one projections onto GP minimizers. 
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A consequence of the Krein-Milman Theorem [IJ vol. 2, Thm. 25.12] is that given any 
7 G r and given any e > there are finitely many GP minimizers (pi and positive coefficients 
Aj (with Xi = 1) such that 



with TrlA^I < e. That is, every element of T can be approximated by a finite convex 
combination of GP minimizers. We also note that part (iii) of Theorem |^ follows from part 
(ii) using Choquet's Theorem |31 vol. 2, Thm. 27.6]. We shall, however, prove part (iii) 
(and Eq. ((TT|)) directly in Section 3 (see Step 4). 

Eq. ()1U() reflects the spontaneous symmetry breaking that occurs in the system under 
consideration. Consider the case of an external potential V{x) which is axially symmetric, 
with symmetry axis given by the angular velocity vector fi. In general, the non-uniqueness 
of the GP minimizer stems from the appearance of quantized vortices, which break the axial 
symmetry, and hence lead to a whole continuum of GP minimizers j22( \2'A\ [21 El IHl 13 ^ • 
Uniqueness of the GP minimizer can be restored by perturbing the one-particle Hamiltonian 
Hq in such a way as to break the symmetry and to favor one of the minimizers, e.g., by 
introducing a slightly asymmetric trap potential V{x). This then leads to complete BEC, 
as can be seen from our Theorem |21 which does not assume any particular symmetry of 
V{x). Note that in the case of a unique GP minimizer. Theorem |21 implies that the reduced 
one-particle density matrix of any approximate ground state converges to the projection 
onto this unique GP minimizer, since Text (and hence F) consists of only one element in this 
case. 

The situation of a dilute rotating Bose gas described in this section contrasts with the 
situation of the absolute ground state of -f/^, i.e., the lowest eigenvalue and corresponding 
state without imposing symmetry restrictions on the wavefunctions. In (23] it was shown 
that Eq. Q does not hold, in general, for the absolute ground state energy. The energy 
per particle in this case is given by minimizing a functional similar to ©, but which now 
depends on one-particle density matrices rather than on wave functions (j){x). In j221 123j it 
was shown that the corresponding energy is strictly lower than E'^^{a) for a large enough 
(and 7^ 0). The density matrix functional has a unique minimizer for any value of 17 and 
a, and in general this minimizer will not be rank one. An analogue of Theorem |21 also holds 
for the absolute ground state. As shown in [221, T consists of only one element in this case, 
namely the unique minimizer of the density matrix functional just mentioned. This implies, 
in particular, that there is no spontaneous symmetry breaking in the absolute ground state. 
We refer the reader to |2S| for more details. 

In the remainder of this paper, we present the proof of Theorems ^ and El 

Acknowledgment. We are grateful to Lev Pitaevskii for drawing our attention to the 
problem of the correctness of the GP equation for a rapidly rotating Bose gas in an email 
correspondence in 1999. 
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2 Proof of Theorem [T] 

STEP 1. Reduction of the Number of Particles to Ensure a Bounded 
Energy per Particle. 

One of the problems we shall face in our analysis is to control three-body collisions, i.e., to 
show that the ground state wave function is suitably small when three particles are close 
together. We have found a way to do this (see Step 4) with the help of a bound on the change 
in energy when three particles are added to the system. It is not evident that this bound is 
always satisfied (although it must be satisfied on average since the total energy is bounded by 
N) and the discussion in this subsection shows how to circumvent this annoyance. If another 
way could be found to control the three-body amplitude or to control the incremental energy 
then the analysis in this section would not be needed. 

Let us consider the Hamiltonian for M < N particles (but still with interaction 
potential depending on A^): 

M 

^M,^ = E^f^+ E ^J^-^j)- (12) 

i=l ^<i<j<M 

This operator acts naturally on all of L'^{E?^'^). We denote the ground state energy in the 
bosonic sector by Eq{M, N). Our goal is a good lower bound on Eq{N, N). 

Let M = M{N) be the largest integer < satisfying two conditions: a.) N — M is 
divisible by 3 and b.) Eo{M,N)-Eo{M-3,N) <6E^^{a). Then E{M +3, N)-E{M , N) > 
6E^^{a), E{M + 6,N) -E{M + 3,N) > 6E^^{a), etc., whence 

Eo{N,N) > Eo{M,N) + 2{N - M)E^^{a) . (13) 

We will prove the following in the remainder of this section. 

Proposition 1. Fix Z > 0, and let Mj and Nj be two sequences of integers, with Mj < Nj, 
limj-^oo Mj = oo and limj^oo Nj = oo, such that Eo^Mj, Nj) — E^^Mj — 3, Nj) < 3Z for all 
j and limj^oo Mj/Nj = A for some < A < 1. Then 

limM ^Eo{Mj,Nj) > XE^'^{Xa) . (14) 

j->oo Nj 

Note that (fT4|) does not depend on Z. 

It is now useful to note that the energy is concave in a (as an infimum over affine 

functions) and thus satisfies 

E^^{Xa) > (1 - A)S^P(O) + XE^^{a) > XE^^{a) . (15) 

The last inequality in (fT^ follows from £^"-^^(0) > 0. 

The sequence M[N) defined above will have a subsequence such that M{Nj)/Nj — > A 
as J — > oo for some < A < 1. If we combine with Z = 2E'^^ {a) we find for this 

sequence that 

lim Eo{Nj,Nj)/Nj > X'^E^^ (a) + 2{1 - X)E^^ (a) = [I + {1 - X)^]E^^ (a) > E^^{a), (16) 

j— >oo 
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which proves Q for this sequence Nj . (Here and in the following, we denote lim inf by lim 
for short.) Together with the upper bound © we also conclude from ()16() that A = 1. That 
is, for Z > 2E'^^{a) the sequence M{N)/N has only 1 as a limit point, and hence ()16() holds 
for the full sequence N = 1, 2, 3 .... 

Our goal in the rest of this section is to prove Proposition which then proves 0, as 
just explained. 

STEP 2. The Generalized Dyson Lemma. 

To get a lower bound on Eq{M, N), we start by deriving a lower bound on the Hamiltonian 
jv ' using Corollary 1 in • This corollary, which is a generalization of Lemma 1 in J2] 
which, in turn, stems from Lemma 1 in Dyson's paper [S], asserts the following. (Note that 
the range of the potential Vj^ is Rq/N, and its scattering length is a/N. We use the "hat" 
to denote Fourier transform.) 

Lemma 1. Let R > Rq/N . Let x{p) a radial function such that < x{p) ^ 1 '^'^^ such 
that h[x) = (1 — x){x) is hounded and integrable (which implies that xip) ^ 1 as \p\ oo). 
Let 

fnix) = sup \h{x -y) - h{x)\ , (17) 
\y\<R 

and 2 

wr{x) = ^fR{x) I fR{y)dy. (18) 

Let Ur{x) be any positive, radial function that vanishes outside the annulus Rq/N < \x\ < R, 
with /]jj3 Uji{x) dx = 47r. Let e > 0. If yi, . . . ,yn denote n fixed points in M^, with \yi — yj\ > 
2R for all i ^ j, then we have the operator inequality on L'^(M^^) 

n n 

-Vx(p)'V + i^7;^(x-y,) > j;((l-e)^C/,j(x-yO-^^K(^-yi)) • (19) 

i=l i=l 

The sums in ()19() are multiplication operators, i.e., they are just functions of x. The 
operator — Vx(p)^V is just the positive multiplication by p^xip)"^ in Fourier space. The 
original Lemma 1 in [T^ has x{p) = 1 ^'^^ h = wr = f = e = 0. 

Clarification: What Lemma^really says is that we can replace the unpleasant interaction 
potential v^^ (which possibly contains an infinite hard core) by a small, smooth, but longer 
ranged potential whose main part, Ur, is positive. There are two prices that have to be 
paid for this luxury. One is to forego a piece of the positive kinetic energy, — Vx(p)^V. The 
second is that the potential is really only a 'nearest neighbor' potential. That is to say, 
the particle at x is allowed to interact with only one other particle at a time. This is seen 
from the requirement that the interaction Ur has range R, but the other particles must be 
separated by a distance 2R. In order to utilize the coherent state inequalities later on in 
Step 3 we have to extend our Ur to an ordinary two-body potential, i.e., we have to be able 
to drop the 2R separation requirement. To do so will require an estimation of the amplitude 
(in the exact, original ground state wave function) of finding three or more particles within 



EHLRS November 22, 2005- GP-rot 



9 



a distance 2R of each other. Clearly, this amphtude is smah, but we find that we have to 
resort to path integrals (or, more precisely, the Trotter product formula) to estimate it. This 
will be done in Step 4 below. 

As an immediate corollary of Lemma ^ we can omit the condition \yi — yj\ > 2R and 
replace (fT^ by 

n 

- Vx(p)'V + i^t;^(x-yi) 

i=l 

n 

^ E [(1 - - vi) - -^^d^ - vi)] U ^ilyk - yi\ - m (20) 

for any set of points yj E M^. Here, 6 denotes the Heaviside step function, given by 6{t) = 1 
if t > and 6{t) = if t < 0. That is to say, if there are only n' < n of the yi that are a 
distance > 2R from all the other y^ then we simply apply (|19j) to these n' coordinates. The 
right side of (|2nj) does not contain the other values of i because the Y\ ^ factor vanishes for 
those. The left side does contain these unwanted yi but, since is non-negative, this does 
no harm to the inequality (|2()|) . 

We apply H2U() to each particle, considering the other M —1 particles as fixed, and obtain 

M 

H,,^^ > (-V.(l - X{p^f)^i + 2pi ■ A{x,) + A{x,f + V{xi)) 
1=1 

M 

+ 2^ [(1 - ^hN-^UR{xi - Xj) - a{Ne)-^WR{xi - Xj)] J] e{\xk -Xj\- 2R) . (21) 

i=l j^i k^i,j 

For the negative part of the interaction (containing wr), we can simply use — ^ fo'^ ^ 
lower bound. For the positive part (containing Ur), we will use the fact that 

J] e{\xk - Xj\ - 2i?) > 1 - ^ e{2R - \xk - Xj\) , (22) 

which follows from the simple inequality Ylj{^ ~ ^j) — ^ ~ Ylj when < Sj < 1 for all j. 

We now use (|21() and (|22() in the following way. We begin by defining a new Af-particle 
Hamiltonian, K, by 

M 

K = Yk^^+ yi 2{l-e)aN-'UR{x,-x,), (23) 

1=1 l<i<j<M 

where Kq is a one-body Hamiltonian to be described next. If Kq were simply (— iV + A)^ + l/ 
then H23|) would be the conventional Hamiltonian with two-body interaction 2Ur. (The factor 
2 arises because each pair i,j appears twice in (PT|) .) 

Unfortunately, Kq has to be a little more complicated because we used up part of the 
kinetic energy in replacing v^^ by Ur via Lemma ^ Pick some r] > 0, and let 

Kq = -V(l - x(p)^) V - 2r]A + 2p ■ A{x) + A{xf + V{x) + r]\x\^ - K{r]) . (24) 
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The constant K(r/) is chosen so that Kq > 0. It is a matter of convenience to include it in 
the definition of K^. It is defined by 

K(r/) = inf spec [— r?A + 2p • j4(2;) + r/lxl"^] . (25) 

The reason for adding the terms — 2r/A and ri\x\'^ to Ko is to ensure that Kq is bounded 
from below and has compact resolvent, and so that K{ri) is finite. (Note: the exponent 4 in 
could be replaced by any exponent > 2 for our purposes. This is due to the fact that 
we have a vector potential A{x) in mind that is bounded by (const. )|x|, as in the case of 
pure rotation. If this is not so (because an external magnetic field has been added) some 
polynomial of higher order than |x|^ could be needed, but our analysis would continue to go 
through.) 

Since there is a 2r/A in p4j) and not just r]A we have that Kq > —r]A + V{x) > —r]A > 0, 
since V{x) > by assumption. This will be convenient later. 

Let {■)i3 denote the M-particle bosonic ground state expectation for the original Hamil- 
tonian Hj^,^ ^. Actually, it is convenient to take the zero temperature limit of the Gibbs state, 
which means that in case of a ground state degeneracy of j^, we would take {■)b to be the 
uniform average over all ground states. Then Eq{M, N) = { H^^j ^ )g and we have, therefore, 
using (EU-lEni), 



Eq{M,N) > inf specif + MK(r/) - r/M(|xi|^)g - 2r/M( - Ai) 



B 



- -j;^{wr{xi - X2))jg —{Ur{xi - X2)9{2R - \x2 - xsl))^ . (26) 

(Note: We made use of the bosonic symmetry to replace Aj by MAi, for example.) 

The term (— Ai)e can be bounded as follows. We have < 2{p + A)'^ + 2A^, and hence, 
using positivity of the interaction potential v^^ , 

M(- Ai>g < 2EQ{M,N) + ^\nfM{\x,f)^. (27) 

To prove Proposition ^ we have to bound the various terms in (|26j) and (|27j) . and that is 
what we do in the following steps. The main term to bound is inf specif, the ground state 
energy of the 'effective Hamiltonian' (|23|) . 

The momentum cutoff x{p) iii dJ will be chosen as follows. Let £{p) be an infinitely 
differentiable, spherically symmetric function with i{p) = for \p\ < 1, i{p) = 1 for \p\ > 2 
and < i{p) < 1 in-between. Then, for some adjustable parameter s to be determined later, 
we choose 

X{p) = i{sp) . (28) 

The potential wr{x) defined in (|1<S|) is then a smooth and rapidly decreasing function that 
depends only on the ratio R/s. It is easy to see that 

/ wji{x)dx < const. (29) 
as long as ^ < const, s. We will, in fact, choose R <^ s. 
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Finally, we are still free to make a choice for the function Ur{x) in Lemma ^ We choose 
it to be a 'hat' function: 

^^(") = \0 otherwise, 

assuming that R > 2^^^Ro/N , a condition that will be amply satisfied by our choice A^~^/^ ^ 
R ^ A^~2/3 ig^^gj. QYL. We remark that the exact form of Ur{x) is unimportant in what is to 
come. We will need only the properties that / UR{x)dx = Att and that ||f/j?||oo < const. R~^ 
for R > Ro/N. 

STEP 3. Coherent State Method for the Ground State. 

We begin our analysis of (|26|) by bounding the main term, inf spec K. This will be done with 
the aid of coherent states, exploiting ideas in [18.^, and is, perhaps, the most methodologically 
novel part of our work. 

The one-body operator Kq has purely discrete spectrum and can be written in terms of 
its eigenvalues ej and orthonormal eigenfunctions \ipj) as Kq = Ylj>i Recall that 

Kq > —r]A + V{x) > —7] A > 0, so ej > 0. We assume the sequence ej to be ordered, i.e., 
Cj+i > ej for all j. For simplicity, we introduce the notation 

W{xi - X2) = {l- e)aN~^UR{xi - X2) ■ (31) 

The well known second quantization formalism involves the operators aj and which 
are the creation and annihilation operators of a particle in the state Ifj). They satisfy the 
usual canonical commutation relations [aj,aj] = Sij, etc. The second quantized version of 
(ESI is 

j>l ijkl 

where Wijki = {^i ipj\W\ipk ^fi)- The operator K. acts on the bosonic Fock space, J^, 
consisting of a direct sum over all particle number sectors. We are interested in a lower 
bound to the ground state energy of /C in the sector of particle number M . Hence we can 
add a term '^j'^j ~ A^)^ to K, without changing this energy. We can then look for a lower 
bound irrespective of particle number. I.e., for any C > 0, we have that inf spec /C for M 
particles is > inf spec /C on the full Fock space, where 

^ = E ^j-^^j + E M'^k^i "^^3^1 + ^ ( E - ' • (33) 

j>l ijkl J>1 

The choice of C will be made later. 

The Fock space T can be thought of as the tensor product of the Fock spaces generated 
by each mode ipy We choose some integer J ^ 1 (to be determined later) and split the Fock 
space into two parts, namely T = T'^ J^^, where T'^ is the tensor product of the Fock 
spaces generated by all the modes i-pj with j < J and where J^-' is that generated by all the 
other modes. 
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Next, we introduce coherent states ^U] for all the modes j < J. (By coherent states 
we mean ordinary canonical Schrodinger, Bargmann, Glauber, coherent states.) The modes 
with j > J will not be omitted, but they will be treated differently from the j < J modes. Let 
z = {zi, . . . , zj) denote a vector in C"^. Let also n(z) denote the projection onto the coherent 
state • -(^zj) G !F^. The symbol jzi®- • -(dzj) is shorthand for \zi)i^\z2)'^- • -(dlzj), and 
\zj) denotes the coherent state for the j*^ mode given by \zj) = e'xp[—\zj\'^ /l+ZjOj] |vacuum). 

The Hamiltonian K. in (|33|) can now be written as 

IC = j dzU{z)(S)U{z), (34) 

where U{z) is an operator acting on J^-^ . The operator U{z) depends on z since it is also an 
upper symbol for the modes j < J. The integration measure is dz = Y[j<j ^^j^Uj with 
Zj = Xj + iyj. This is discussed in j^^j. As an example, the upper symbol for Oj is Zj 
and for it is zj, but for a^a^- it is — 1. Thus, to a term such as alojai^ai with i,j < J 
and k,l > J would correspond the upper symbol operator Zi Zj a^a^ . 

It is easier to compute the lower symbol (which is denoted by u{z)) than the upper 
symbol U{z). It is obtained simply by replacing a)- by Zj and by Zj in all (normal- 
ordered) polynomials, even in higher polynomials such as a)-a- or o)-a}-a-a- . An equivalent 
definition of the lower symbol of any polynomial V in the a^ 's and aj's (normal-ordered or 
not) is the expectation value u{z) = (zi (g) • • • (gi zj\'P\zi ® ■ ■ ■ ® zj). In the case considered 
here, u(z) = [zi® ■ ■ ■ ® z j\K,\zi ® ■ ■ ■ ® zj) . 

The lower symbol is useful because the upper symbol can conveniently be obtained from 
it as dn] 

U{z) = e-^-^-u{z) = u{z) - d^diu{z) + i [d^dif u{z) , (35) 

where d^dz = Y2j<j (^zjdzy (In the general case there would be higher order derivatives on 
the right side of (|35|) . but not in our case since m(z) is a polynomial of order four.) Note 
that (|34|) implies that 

inf spec /C > inf (inf spec C/(z)) , (36) 

z 

since / dzn(z) = Ijr< and n(z) (g) C/(z) > [inf spec [/(z)] n(z). 

Our goal in the rest of this subsection is to derive a lower bound to inf spec C/(z) for a 
fixed z. The reader might wonder why we use coherent states only for modes j < J and not 
for all modes. The reason is that the upper symbol for the operator eja)ja- is e^d^jp — 1), 
and the —1 term is a term that we do not want when minimizing for a fixed z. We make an 
error in the energy of the form — X^j<jej and for this reason we cannot take J = oo. But 
we can, and will let J ^ oo as ^ oo. 

3a. Lower Bound on the Lower Symbol n(z). 

In order to derive a lower bound to U (z) and the bottom of its spectrum, we start by deriving 
a lower bound to the lower symbol u{z), which is the first term in H35|) . This symbol can 
be conveniently expressed in terms of the function $z £ L^(R^), parametrized by z G C'^, 
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given by 



Note that ||^>z||2 



i<j<J 



(37) 



j<j 



Denoting T = Y.k>j '^kolaf^, we have 



Zl 



zj 



(38) 



There is a mild abuse of notation here, which will continue for the rest of this paper, and 
which we hope will not cause any confusion. The operator • e,ato ■ acts on T while the 
vector I Zl (gi • • • (8i zj) is in , so the left side of (|5S|) defines an operator on J^-^ in an obvious 
way (actually, it defines a quadratic form). The right side must also be an operator on J^'^ , 
and it is so if the number {^z\Ko\^z) is regarded as a number times the identity on J^-^ . 
Similarly, with A^^ = denoting the number of particles in the modes > J, 



Zl 



O ZJ 



fiV> + ||$z||2-M)^ + ||$,"2 



-z|l2 



> (ll^zlli -M)^ -2e7^Mr. (39) 



Here, we used the normal ordering Ej<jaj«j]^ = Y.i<jY.j<j + Ylj<j'^](^j^ Al- 

lowed by the elementary bound < T jej. 

The interaction part of n(z) is obtained by replacing by Zj and by Zj when j < J. 
We will now derive a lower bound on this term. It is convenient to introduce the notation 



(40) 



Since > 0, it is possible to neglect the interaction between modes > J for a lower bound. 
More precisely, let P = Yli<i<j Q = 1 — P. The two-body operator W{x — y) 

is then bounded from below by 



W = {iP + Q)<S)iP + Q))WiiP + Q) (P + Q)) 

> (P 8) P)W{P (g)P) + {P^ P)W {P®Q + Q®P + Q®Q) 
+ iP0Q + Q<S)P + Q<S)Q) WiP <S) P) , 



(41) 



since the missing term on the right side of (|1T|) is {Q <Ei Q + P Q + Q <^ P)W{Q <^ Q + P <^ 
Q + Q <^ P) >0. We thus have that 



Zl 



Zl 



(gl ZJ 



I 



> I(^>z) ^z\W\ipk g) 'Pi)a^ai + {ifk ® ipi\W\^z ® ^^)a\cL 

>J kl>J 

2 Y^^^ ® ^zlV^l^z (S) iPk)af, + 2 ^($z (g ^k\W\^z (g ^z)a 



kl>,J 



(42) 



k>J 



k>J 
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Here we used that W is symmetric, implying that in the last line we could replace |<I>z 
^Pk + ^kf^ ^z) by 2|$z (g> ipk). 

We seek a lower bound to the last two expressions in H42|) . Note that, for a general 
operator A, \A + = + A^"^ + AA^^ + A"!"^ < 2A"I'A + 2AA^ by Schwarz's inequality, 
and so 

{A + A^f<A\A\'^ + 2[A,A^]. (43) 

We apply this to the second line in with A = '}2ki>j '^kio-j^a^ and = {(^z®^7,\W\(pk® 
ipi). The commutator is 

[A,A^] =4 Yl (^Ul^^'lrn + 2 J] |Cfc,P . (44) 
klm>J kl>J 

The last term in ()44() is bounded by 

5^ Icfcd' < I^'^'l' = ® ^zlVF'l^z ® ^-z) < ||T^||oo/(^z) . (45) 
kl>J kl>l 



The first term on the right side of (|44() can be bounded as 

4 



XI 0.1(^1 CkmClra < Y ^ "'k^l'^''"^'^^"' - ^ II W^lll II V^>z II2 ^ • (46) 



klm>J m>lkl>J 



This can be seen as follows. The integral kernel a of the one-particle operator defined by 
the matrix Y^^^^q^cim is given by 



a{x, x') = J dy\^,{y)\^Wix - y)Wix' - y)<^.{x)<^,{x') . (47) 
Using Young's and Schwarz's inequalities, we have, for any function / on M'^, 



dxdx'f{x)f{x')a{x,x') < J dx dx' dy\^^{y)\^W{x - y)W{x' - y)\<^^{x)f{x)f 

< ||W||?||fz||^||<l>z/||i < ||W^||?||^z||^||/||^ (48) 

Hence (jMl follows by applying the Sobolev inequality ||/||| < (4/3)(27r2)"2/3||V/||| both to 
/ and to <I>Z5 and using the fact that — A < ti'^Kq. 

To get an upper bound on \A\'^ we use Schwarz's inequality again to obtain 

1^1^ - \Y1 ~~ ) ( ^ ^rnenai^aia^aA (49) 

\kl>J ^ ' / \mn>J J 

for any sequence of positive numbers Cj. We choose Cj to be the eigenvalues of Kq, in which 
case ^ 

mn>J \k>J / 

Moreover, 



^z ® ^z ) . (51) 
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We have the following two operator inequalities; the first comes from the fact that Kq > ej 
on the range of the projector Q and the second comes from Kq > —rjA: 

Q 2 2 

— < < . (52) 

Kq - Ko + ej - -r]A + ej 

Denoting the integral kernel of (—A + fi)^^ by 

1 g-v^|a;-a;'| 

^.(--^') = ^^3^, (53) 
we see that (|^T|) is bounded above by 



dxdydx'dy' ^^{x)^^{y)W(x - y)Kj/^{x - x')k^j/^{y - y')^^{x')^^{y')W{x' - y') 
dxdydx'dy' |$z(x)P|<i>z(2/)PVF(x - y)K^/r,{x - x')Kj/^{y - y')W{x' - y') 



4 

< - 

Here, we used Young's inequality for the {x',y') integration, as well as the fact that \\k^\\2 = 
By putting all this together, we have that 

{A + At)2 < 4||T^|U/($z) + ^V-'\\W\\l\\V^.\\l T + _i^&/(ci>,) T' . (55) 

Since the square root preserves operator monotonicity, we can take the square root on both 
sides of (|55p. By the triangle inequality, we can take the sum of the square roots of each 
term on the right side. Finally, applying the Schwarz inequality to the first and third term, 
we conclude that, for any 6 > 0, 

\A + At| < <5/($,) + ^\\w\\oo + ^yj??~'/'l|VF||i||v$z||iVr 

+ ej'/' (\\WhI{<f,) + (2^/2)-^) T . (56) 



We now proceed similarly with the last term in (|42|) which is linear in and a|, . Denoting 
Cfc = (^z i8) $z|W^|^z ® V'fc); we have that 

\k>J J \k>J ^ / \k>J ) k>J 

Using Holder's and Sobolev's inequality, 

J^|cfc|2= f dxdydz\^,{x)\^\<^,{y)\^\<!>,{z)\^W{x - y)W{x - z) 
k>l 

<l|W^ll3/2ll^z||^/(<^z) < ||V$z||i/($z) . (58) 
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Moreover, using (|52() again, together with Young's and Sobolev's inequahties, as well as the 
fact that the 3/2 norm of /c^ is given by 2^^^^ fi~^/'^ /3, we find that 

^ kfcP < 2 /• ^^^y^^/^y/^;(^|^^(^/)|2^(^ _ x')k,^/^{x - y)My)\^,{y')\^W{y - y') 



k>J 

2 



<-j dxdydx'dy' |$,(x)n$,(x')ri^(x - x')K^,r,{x - y)\^.{yTW{y - y') 



This implies that 



\k>J J 



S 5p;^l»»'ll^''ll»'ll?"llV*ji/(*.) + ^,-"=l^l|V<|.jP(*.)T (60) 



< 



2 

3 (27^2)1/3 II "''11- " "1 ^ 3^2/3 



again using the triangle and the Schwarz inequality. As mentioned above, operator mono- 
tonicity is preserved by the square root, and hence we can take the square root on both sides 
ofEq. (|Hn|). 

This completes the lower bound on the lower symbol ti(z). For the convenience of the 
reader, we repeat the bound just derived: 

n(z) >{^^\Kq\^.) + /(CD,) {l-5- e}'/^||Ty||ir) + - Mf (61) 

|2 ^ T(S>\\ i 9, /? I IIWI|1/6||H/||V3 , ^„-l/4||w||l/2 -1/4 . ^ 



We note that in the following we will choose J large enough so that the last term in 1)61^ is 
positive and thus can be neglected for a lower bound. (Recall that ej — > oo as J — > oo.) 

3b. Lower Bound on the Remaining Terms in U{z). 

A lower bound on the first term on the right side of psp is given in ()6ip and, therefore, to 
get a lower bound on the upper symbol U{z), it remains to bound the last two terms on 
the right side of H35|) . The very last term is positive, as will be shown now, and can thus be 
neglected for a lower bound. Namely, 

^(d.d^fuiz) = lid,d^f ((^z®^'z|v^^|^z®^z) + ^||$z||2) (62) 

c 

= ^ X] iv'i'^ "^j + y^i\w\(pi'^ (Pj + ifj 1^ ifi) + —J{J + 1) >o . 

^<ij<J 
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The remaining expression, dzdz,u{z), consists of the following terms. First, from the 
one-body part (1^5]) of the Hamiltonian we obtain a contribution Ylj<j^j- Second, from the 
term H39|) (see also (jSSl) that was introduced in order to control the particle number, we get 

- {{2N> -2M + 1 + 2\\^,f)J + 2||cl>,f} < _(J + 1)||<I>,||2 + _ (2iV> + i) . (63) 

Finally, the following three contributions are obtained from the interaction part. From 
the part where all four indices are < J, we have 

^($z ® (fj + ifj $z|W^|$z ® + fj ® ^z> < 4^($z ® V^ilW^l^z ® ^j) 

<4j^||W||i||$,||2||^,||2 < (4/3)3/2 -i^r?-i/2||VF||i||VcI>,||2j^^. (64) 

Here, we used the inequalities of Young, Holder and Sobolev as well as the facts that —A < 
r]~^KQ and {ipj\KQ\ipj) = ej in the last step. From the term with 3 indices < J, we get 

2 ^ ^($z + ® ^zlM^lv'i O + adjoint . (65) 

j<Jk>J 

Using ((S7|l. this time with = 1, is bounded above, as an operator, by 

1/2 



k>J 



(66) 



Similarly to (|58|) . we can derive the bound 

^ fk\W\<!>^ ^^i + ifi^ $z)|^ < 4||Ty||3/2l|W^I|l||^z||6ll'/'i|l6ll'/'i||2 • (67) 

k>l 

Since \\ipi\\2 = 1 and < (4/3)(27r2)-2/3||Vy,^||2 < (4/3)(27r2)-2/3^-ig.^ this implies that 

® - 8(V3)^/^^^^||T^||f ||Ty||^^||VcI>z|br?-^/^5:ef ^iV> + 1 



i<J 



- 4(V3)^/^^^^l|T^llf llM^II^/^-^/^E^^^^^ (l|V$z||i + iV> + 1) . (68) 

Here, Schwarz's inequality was used in the last step. 

The last term to estimate is the one coming from 2 indices < J, given by 

{y^j ^y^k + Vk® y^j\W\ipj ®Lpi + ipi® (Pj)alai 

j<Jk,l>J 

< (4/3)3/2^r?-3/2||H^||i^^r. (69) 
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This inequality can be seen as follows. For any one-particle function /, 

{^i(^f + f(S)^i\W\^i(S)f + f(E)ipi)<Mipi(E)f\W\^i(S)f) < 4||PF||i||/||^||(/p,||i 

< (4/3)3/2 -L^-i/2||^||^||v/||2^. (70) 

The last inequality is the same is in H64() . The result now follows using — A < r]~^KQ. 
Altogether, we have thus shown that 

d^d-M^) < Y^e. + ^(j + im.wi + ^-§f{N> + 1) 

i<J 

+ mf'l^.V-'^'WWh iim.Wl + V-'T) (71) 



+ 4(4/3)^/^ ^^-i^||l^||f||T^||VV^/^ E ^f" + N> + 1) 



i<j 



This finishes our lower bound on the upper symbol U{z). To summarize, we have shown the 
following operator lower bound to the operator U{z): 

U{z) > right side of dHU - right side of (HH) . (72) 



3c. c-Number Bound on T. 

We are interested in the ground state energy of U{z) for a fixed z G C"'. Since T and 
are the only operators appearing in (|H1|) and (|71|). this quantity can be bounded from below 
using (|H1|) and (|71|) if we can evaluate the expectation values of T and N'' in the ground 
state (or one of the ground states) of U{z). 

Let ( • )z denote the expectation value in a ground state of U{z). We can use two simple 
facts: i.) Since \/T enters (jblf) negatively, we can use the concavity of the square root to 
replace {VT)z by ^/'(T)^ for a lower bound, ii.) Since appears positively in H71() . and 
hence negatively in (|72|) . we can replace it by the upper bound < T/ej. 

For the purpose of bounding {T)z we can use a lower bound to U{z) that is much simpler 
than H72|l . This is obtained by totally neglecting both the interaction part and the part 
controlling the particle number in u{z). These give positive contributions to u(z) (since 
u{z) is the expectation value of fC in the coherent state). We have to be more careful 
about dzdzU^z), however, because this contains some negative terms, as given in ((7T|) . (The 
annoying fact is that an upper symbol of a positive operator need not be positive, although 
the lower symbol is always positive.) 

Proceeding in the manner just described we have that 

C/(z) > ($z|i^o|$z) +T- d^dM^) . (73) 

Now let us estimate the various terms in dzdzU^z) in (fTT]) . We have ??||V^>z||2 < (^z|-f^o|*l'z)- 
Also ll^zlli < ($z|i^o|^'z)/inf spec(-7?A + y(x)). Moreover, \\W\\i < Ana/N, and ||T^||oo < 
6a/{R^N). 
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We will choose R ^ A^^^/^ below. Therefore, ||Vl^||oo ^ The operator N-^ can be 
bounded in terms of T as < T/ej. Note also that M = 0{N) by assumption. Hence we 
see from (f7T|l and (f75)) that, for N large enough (depending on the parameters rj, C and J), 

U{z) >\T- const. , (74) 

where the constant depends only on r/, C and J, but not on M or A?". 

The value of ()74() is that it allows us to control the value of {T)^^, and thereby control 
infzinf specC/(z), which is our lower bound to the ground state energy of tC. There is some 
number independent of all parameters, such that infzinf spec [/(z) < ME/2 — const, 
because infz inf spec f/(z) is less than the known upper bound to the ground state energy of 
/C. Then we can, and will restrict our attention to z's with (T)z < ME because only those 
values of z are relevant for computing inf ^ inf spec C/(z), as (|74|) shows. Only the existence 
of E and not its value is important. 

We conclude from ()72() and the fact that {T)^ < ME for the z in question that 

inf inf spec C/(z) > inf£:[^>] , (75) 

z ^ 



where 

Sm = imom + Di /(CD) + ^ (||$||2 - Mf - D2\\Vm -D^- ^(J + . (76) 

The notation is the following: 



4 



37r2/3 

+ ^Mv-'^'\\WhM'/'E'/' + (4/3)3/2^^-V2||^||^ ^ ^ 

+ 4(4/3)^/^ ^^-i^inifinivv^/^ , (78) 

and 

^3 =E^^ + ^ (^^/^^ + ^) + (4/3)=^/^^^-^/^||W^||iMi? 

i<J i<J 

+ 4(4/3)^/^ ^^-i^llT^llf IIT^IIVV^'/^ (^^^:^' + ^) + Jll^lloo • (79) 



We have neglected the last term in (|61|) containing (1 — ej^''^(27r^r/^/2) i — 2C/ej), assuming 
J to be large enough to make this term positive. (Recall that ej — > oo as J — > cxd.) 
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Our final result in this section, ((75|) - H76|) . might not appear to be useful at first sight, 
but the reader should note that the first two terms in H76|) are essentially the GP energy 
expression. The term ($|i('o|^) is the relevant (i.e., low momentum) part of the kinetic 
energy f \{i\7 — A)^\'^. The coefficient Di equals 1 to leading order and /($) is essentially 
the GP quartic term Ana / l^l"^ (up to errors which will be controlled). Moreover, for C 
large enough the term C(||$||| — M)'^ /M ensures that we have the right particle number. 
For an appropriate choice of the parameters J, rj and R all other terms are of lower order 
as — > oo, as we shall show. 

STEP 4. Bounds on Three-Particle Density. 

So far we have bounded the main term in H26|). namely inf specif. Of the various other 
terms in ()26() that have to be bounded, the one that is most intuitively negligible, but which 
we find the hardest to control is the last term in H26|). To show that it is small we have to 
show that the probability of finding three particles within a distance 2R of each other (in a 
true ground state of H^^j ^) is small. This is accomplished in this section. 

We begin with a lemma about the possible size of the expectation value of a function of 
the coordinates of three bosons. Recall from Step 2 that { ■ )b denotes expectation value in 
the bosonic, zero-temperature state of the M-body Hamiltonian H^j 

Lemma 2. Let ^{xi,X2,X'i) be any positive function of xi, X2 and X3 G M^. With V = 
the one-body potential appearing in H^^j^, we define the three-body, independent particle 
Hamiltonian 

= -Ai - A2 - A3 + V{xi) + V{X2) + F(X3) . 

Let a > and let e~"'^{xi, X2, X3 ; yi, 2/21 2/3) be the 'heat kerneV of h at 'inverse temperature' 
a. Finally, consider the modified integral kernel 



e 



-ah I 



xi,x2,x3 ; yi, y2, y3)V^ixi^x2,x3)Ciyi,y2, ya) (80) 



and let A denote its largest eigenvalue (i.e., its norm as a map from L (M.^) to L (R-^) ). Then 

(e(xi, X2, X3)>g < Aexp{a(So(M, N) - Eo{M - 3, N))} . (81) 

Note that for the M and N under consideration here, we have Eq{M, N)—Eq{M—2>, N) < 
3Z, as explained in Step 1. It is the appearance of the peculiar difference Eq{M, N)—Eq{M— 
3, N) in Lemma[2lthat led us to the discussion in Step 1. If the three-body correlations could 
be bounded more expeditiously than is done here. Step 1 could be simplified. 

Proof. We denote by Tr [•] the trace over all of L^(R^*^), not just the bosonic states, and 
by Pb the projection onto the bosonic (i.e., symmetric) subspace. Note that exp{— ^} 
is trace class for large enough /?, by our assumption on the logarithmic increase of the 
potential V{x). (This follows from the Feynman-Kac-Ito formula, together with the results 
in the appendix.) Hence 

Tr r^p"""^M,jv R 1 
n^oo Tr[e~ «'^Pb] 
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independently of a, of course. Note that Hj^j ^ commutes with Pb so e 



-anH, 



Pb is self- 



adjoint and positive. The multiplication operator ^ is also positive and we can write 
^^-anHj^j j^ _ j^g-°^A/,jvp^jg-°("-i)^Af,]vpj^, Holder's inequality for traces of positive op- 
erators states that TiAB < {TrA"}i/"{TrS"/(""y}("-i)/", and therefore 



(83) 



Consider the bigger projection Pb, which symmetrizes only among particles 4, 5, . . . ,M. 
It commutes with -H^,^^ and also with ^, and hence e~'^^^'-'^ P\, < e~"^^^''^Pb- Since ^ > 0, 
this yields the upper bound 



Tr[e"""-^M.^Pb] 



(84) 



We now claim that 



(85) 



where TVa and Trj\/-3 denote the trace over the first 3 and last M — 3 particles, respectively. 
Taking the limit n ^ oo this proves (|81() . 



To show , we write 



M.N 



f + W, with W denoting the inter- 
action between the first 3 and the last M — 3 particles. Note that Vl^ > 0. Using the Trotter 
product formula, we first replace each factor e""^«,]v by (e-"^3,iv/mg-a(/fM_3,iv+VK)/m^m f^^, 

some integer m. (Here we abuse the notation slightly, omitting to write tensor products and 
identity operators.) For x = {xi,X2,X3), let fe(x, x') denote the integral kernel of e~"^3.^/™'. 
Denoting by TVx the multiplication operator on the subspace of the last M — 3 particles 
obtained by fixing the first 3 to have positions x, and introducing nm integration variables 
Xjj, 1 < i < n, 1 < J < m, we can write 



Tr 



-aH-i^N /m^-a{HM-3,N+W)/m 



-a(/i'M-3,JV+VKx,, )/m 



, (86) 



where we identify Xj(„_,_x) = ^(i+i)i ^n(m.+i) = ^i,!- By Holder's inequality for traces, 
we can estimate 



IVm- 



< sup TrM-3 

ij 



< Tr 



M-3 



Phe 



(87) 



where in the last inequality we used the fact that 14^x ^ and that the partition function 
is monotone in the potential. By the Feynman-Kac-Ito formula 24, Sect. 15], the integral 
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kernel A;(x, x') is bounded in absolute value by the kernel of e "''/™. Using this estimate 
and rewriting the integrals as a trace we obtain 



Tr 



-aH3^N/m,^-a(HM~3,N+W)/m 



Letting m ^ oo this yields 



□ 



We now use Lemma|21to obtain a bound on the various terms in (|26j) and (|27j) . Lemma[21 
immediately implies that 



(89) 



with II • lloo denoting operator norm. For positive operators, the operator norm is bounded 
by the trace, in this case given by Tr\x\'^e°'^'^~^^^^h This expression, in turn, is bounded for 
a large enough, as shown in the appendix. In exactly the same way we can bound (|xi|^)g. 
Moreover, we have that 



{wr{xi - X2))_g < e' 



='"^'|^/i^e--'^V^||oo<e3°^ 



(47ra)3/2 7jg3 



'Wji{x)dx . (90) 



The last inequality can be seen as follows. Denote by k(x,x') the kernel of e"*-^"^^^)^. The 
Feynman-Kac formula implies that k{x,x') < (47ra)^^/^ for any positive V{x). Hence, for 
any function f £ L^' 



dxdx'dydy' f{x, y)\/wR{x - y)k{x, x')k{y, y')\/wB{x' - y')f{x, y) 
< (47rQ)^^/^ / dxdx k{x,x ) 



dyy/wR{x - y)\f{x,y)\ 



< (47ra)"^/2 ( r^R] / dxdx'k{x,x') 



dy\fix,y)\' 



dy^/wR{x' - y)\f{x',y)\ 

1/2 p |. -| 1/2 

dy\f{x',y)[^ 



(91) 



where we used Schwarz's inequality in the last step. The result now follows from the fact 
that e"(^-^(^)) < I. 

Similarly, repeating the above argument with X2 in place of x and (xi, 2:3) in place of y, 
we obtain 



[Ur{xi - X2)e{2R - \X2 - X3|)>„ < e 



3oZ 



1 



UR{x)dx / 9{2R- \x\)dx 



Sir 



(92) 



This finishes our bounds on the various terms appearing in (|26() and (|27|) . 



STEP 5. Collection of All the Terms and the Final Inequality. 

In this section we concatenate the various pieces of the lower bound to the energy Eq{M, N) 
in H26() . and finish the proof of Proposition ^ Inequality (|26|) contains several terms. All 
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except inf spec K were bounded in Step 4 and in (P7)l . The essence of Step 3 is the bound 
on the main term 

inf spec K > inf inf spec U{z) > inf £[^] , (93) 

z <I> 

where S[^] is defined in H76() . 

Let us begin by disposing of the terms mentioned in Step 4. As shown there, < 
const, and < const, for some constant depending only on Z. (Recah that Z is a 

fixed number of order 1.) Moreover, from (|9U|) and H29|) we see that (recalhng that M < N) 

wr[xi - X2))y, < const. (94) 



This term wih thus be neghgible, if i? ^ as ^ oo (keeping e and s fixed for the moment). 
We are free to choose the dependence of R on N, and we choose R to satisfy 

N-^/^ > i? > AT ^ oo . (95) 

The last term to estimate is then 

-j;^{Ur{xi - X2)e{2R - \x2 - xgl))^ < const. aNR^ <C 1 . (96) 

Hence it follows from (|2(i)) and (|27)) that, for any fixed s, £ and rj (recalling that A = 

limjv^oo M/iV), 

lim ^{l + 4r])Eo{M,N) > lim ^ inf + Ak(?7) - const. Ary . (97) 

TV^oo A Ar_^oo A * 

The only thing left is the minimization of £ [^] given in ()76() , which contains the numbers 
Di, D2 and in (f77j) - (f7^ . To evaluate them as A^ ^ 00 we note that < Aira/N, 

and 1 1 1 1 00 ^ A^ for our choice of R in (|^^ . Hence, we see that 

lim lim lim Di = 1 , lim lim D2 = , and lim — D3 = 0. (98) 

Using the fact that both ||V^'||2 and ||<5||2 are bounded relative to (^li^ol^)) and rescaling 
we obtain 

lim lim ^mf£[^] > (99) 

lim inf I A($|A:o|^) + (1 - e)aA2 /" |$(x)|^|$(y)pt//j(x - y)dxdy + CX - iX 

R^o * I J 

Note that the infimum can obviously be restricted to a set of bounded (^>|-fCo|<I?), in- 
dependent of R, since Ur > 0. Since Kq > —rjA this implies that we can assume that 
||V^>||2 is bounded independent of R, and hence also by Sobolev's inequality. Using 

the inequality (proved below) 



J Mx)fMy)fUR{x - y)dxdy - 4^||^>||^ 



< 87ri2||^>||^||V«>||2, (100) 
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we see that we can interchange the hmit and the infimum and thus obtain 

lim lim ^ inf £:[$]> inf |a($|Ko|^) + (1 - e)47raA^||^>||l + CA (||^>||2 - l)^ . (101) 

Ineq. (jlUOj) can be obtained in the following way. Using Schwarz's inequality, as well as 
/ UR{y)dy = Air, 

J mx)\^my)\^UR{x - y)dxdy - 47r||$||l 
< j dyUniy) j dx\^x)\W^{x)\ + \^{x + y)\)\\^{x)\ - \^{x + y)\\ 

<2ml j dyURiy)Q' dx Mx)\ - \<l>ix + y)\\^'^ ' . (102) 

The result now follows from the fact that |||<I>| — |<1>( • +y)|||2 < |?/|||V$||2, which can be seen 
by evaluating the norm in Fourier space, using |1 — e~^v-y\^ < and also using the 

fact that / \jR{y)\y\dy <Rf UR{y)dy = AttR. 
Now, letting C ^ oo, we infer from (|101j) that 

lim lim lim— inff[$]> inf {A($|Ko|^) + (1 - e)4vraA^||$||l| . (103) 

C-^oo J~>oo N^oo N * 11*112=1 

The final step is to remove the momentum cutoff in Kq, i.e., to let s ^ in Eq. (|103|) . 
Again, we claim that we can interchange the limit and the infimum, at least to obtain a 
lower bound. Let $s denote a minimizer of the functional on the right side of (|103p . Since 
Kq > —r]A + V{x) and V{x) ^ oo as |x| ^ oo, a sequence with Sj — > as j — > oo 
lies in a compact subset of L^(R^), and hence there exists a subsequence which converges 
strongly and pointwise almost everywhere (both in p-space and x-space) to a function 
$0 as j ^ oo, with 11*^*0 II 2 = 1- All the s- independent terms in the functional on the right 
side of H103|) are weakly lower semicontinous. Moreover, by Fatou's Lemma |11| . 

lim / /(1-Xs(p)')|$s(p)|'dp> / p^\$oip)\^dp. (104) 

Hence the infimum and the limit s ^ can be interchanged for a lower bound. In combina- 
tion with inequalities (|103|) and (|97|) . we find that 

hm ^{l + 4r])Eo{M,N) 

> inf {\{^\-A + 2p- A{x)+A{xf + V{x)\^) + {1-6)4^0X^11^1} -const. Xr]. 

||$||2=1 

(105) 

(For a lower bound we simply dropped the positive terms —2r]A and t7|x|^ in Kq.) By letting 
r] ^ and e ^ Proposition ^ is proved. As explained in Step 1, this proves Theorem ^ 

Remark about the optimal choice of the parameters: In Eq. (|95() we showed how the 
parameter R has to depend on N, as N ^ oo, in order to obtain the correct limit for the 
energy. The explicit dependence on N of the other parameters J, C, s, r] and e need not 
be specified so closely (unless we wish to obtain a detailed error estimate). It suffices to let 
J oo, C — > cxD, s — > 0, r/ ^ and e — > (in this order) after taking the N ^ oo limit. 
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3 Proof of Theorem [2] 
STEP 1. Proof of Part (i). 

The fact that F is a convex set follows easily from its definition. Namely, if and \ are 
two approximate ground state sequences, and < A < 1, then Xj^^ + (1 — X)% is certainly 
also an approximate ground state sequence, whose reduced one particle density matrix is 
givenby A7« + (1-A)7^i). 

Compactness of T is also not difficult to see. Given a sequence 7^ G F, the Banach- 
Alaoglu Theorem implies the existence of a subsequence such that 7^ ^ 700 for some 700 
in the weak-* sense as z — > cxd. As already remarked in the introduction, the fact that 
Ti Ho'ji < const, implies that 7^ — > 700 in trace norm. To prove compactness we have to 
show that 7oo G F. 

By definition, corresponding to every ji there is an approximate ground state sequence 
7jy -. That is, there is a number Ni such that N > Ni implies that W'jj^^ — < 1/i and 
\N^^Tr Hj^jj^ ^ — £^"^^(0)1 < (Here, || • || denotes trace norm.) We can assume that 

A'j ^ 00 as z — > 00. Now, for given N, let i{N) be the largest integer i such that N > N^. 
Then i{N) ^ 00 as ^ oo, and hence the sequence J^^f^^-^ is an approximate ground state 
sequence. Moreover, 117^^^^ - 7oo|| < - %(jv) II + H(n) - 7oo || ^ as A^ ^ 00. This 

proves that 700 G F, and hence F is compact. 



STEP 2. An Extension of Theorem [H 

A key step in the proof of Theorem |21 is an extension of the lower bound in Theorem ^ 
to the case of a perturbed Hamiltonian, where we replace the one-particle part Hq of the 
Hamiltonian ^ by Hq + S, where S is a bounded hermitian operator on the one-particle 
space L^(M^). Let Hj;^^^ denote the perturbed A-particle operator 

N 

1=1 

and let E^q^\n) = inf spec Hj^^^ denote its ground state energy. Correspondingly, define the 
perturbed CP functional as in ©, with Hq + S in place of Hq, and let E'^^{a) denote 
its infimum over all 4> with \\(j)\\2 = 1- Then we have the following extension of Theorem 
to whose proof we will devote the remainder of this subsection. 

Proposition 2. For all bounded hermitian operators S, 




We start by noting that in order to prove Proposition |21 it suffices to prove it in the 
special case in which 5 is a finite rank operator with exponentially decaying eigenfunctions. 
In particular, we can assume that its integral kernel S{x, y) satisfies a bound 



\S{x,y)\<Bexp{-D{\x\ + \y\)) 



(107) 
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for some positive constants B and D. This can be seen as follows. Let {fi}fZi be an 
orthonormal basis for L^(M^) such that |/i(x)| < Bi exp{—Di\x\) for some choice of constants 
Bi, Di > and let Pn denote the projection onto the first n of these functions. Clearly, 
Pn — > I strongly as n ^ oo. Then, for any bounded S, PnSPn is of the desired form, i.e., it 
has finite rank and its integral kernel satisfies a bound of the form ()1U7() . For any one-particle 
density matrix 7, 



Tx[^{S - PnSPn)] 



< 



^ -.{S- PnSPn)- ^ 



fij \ " rrr 



Tr[ffo7] , (108) 



— 1/2 

with II • II denoting operator norm. Since is compact and, therefore, is the norm limit 

of finite rank operators, it is easy to see that the norm in H1U8() goes to zero as n — > co. On the 
other hand the set of numbers Tr[ffo7] that arise from those 7's that come from approximate 
ground states is bounded. Consequently, both sides of (|1U6() can be approximated to within 
any desired e by replacing S by PnSPn and choosing n large enough — which implies the 
statement. 

Thus we can assume (|107p henceforth. The proof of Proposition ^ then follows exactly 
the same lines as the proof of Theorem ^ In fact, our proof of Theorem ^ has the advantage 
of being almost completely independent of the exact form of the Hamiltonian. The only 
place where we used the explicit form is Lemma [21 which was used to bound expectation 
values of certain one-, two- and three-body operators in the zero-temperature state of H^^^ ^. 
We now have to bound the expectation value of these operators in the zero-temperature 
state of H^jIj, which we denote as ( • (Here, the operator H^f]^ is defined in the obvious 

way. Its ground state energy will be denoted by Eq (M, N).) To this end, Lemma|21can be 
extended in the following way. 

Lemma 3. Let C(xi, X2, 2^3) he any positive function of xi, X2 and X3 G M^. Let S denote 
the rank one operator on the one-particle space with integral kernel given by the right side 
of jlO'^ - With V = the one-body potential appearing in H^,j^, we define the three-body, 
independent particle Hamiltonian 

h^S) ^ _ A2 - A3 + V{x{) + V{X2) + V{x^) -S1-S2-S3. (109) 

Let a > and let e~°'^*^' (xi, X2, 2:3 ; yi,y2,y3) be the 'heat kernel' of /i^"^-* at 'inverse 
temperature' a. Finally, consider the modified integral kernel 



e-" 



''^^\xi,x2,xs ; yi,y2,y3)V^(.^i^^2,x3)^{yi,y2,y3) (110) 



and let A^^^ denote its largest eigenvalue (i.e., its norm as a map from L'^(MP) to L'^{MP)). 
Then 

<e(a;i,a;2,X3)>Jf^ < A(^)exp{a(i?f)(M,iV)-4^)(M-3,iV))}. (Ill) 

The proof follows along the same lines as the proof of Lemma |2l except for one step. 
Before Eq. (|88|) . it was necessary to get an upper bound on the absolute value of the integral 
kernel of exp{—aH,^j^/m} in terms of the kernel of exp{— a/i/m}, which can be obtained 
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with the help of the Feynman-Kac-Ito formula. In the case considered here, we need an 
upper bound on the integral kernel of ex.]){—aHj^/m}. We will now show that the absolute 
value of this kernel is bounded above by the kernel of exp{—ah^^^/m} for the modified 
three-particle operator 

This claim follows from the Trotter product formula, together with the Feynman-Kac-Ito 
formula, in the following way. Since 5* is a bounded (in fact, finite rank) operator, we can 
write 



/m 



lim 



«5 
n 



n/m 



By the Feynman-Kac-Ito formula, (|1U7|) and the definition of 5, 



n ■ 



n/m 



< 



n ■ 



n/m 



ih(^'> /r. 



fll2) 



(113) 



This 



In the limit n — > 00, the operator on the right side converges strongly to e" 
proves our claim. 

For the application of this Lemma, as in Section 2, Step 4, it is necessary to have some 
bounds on the kernel of e~"'^''^\ In particular, we need that the kernel is bounded, and 
that its diagonal decays for large at least like " for some positive constant. As 

for the case 5 = 0, these properties are again shown in the appendix. It is there that the 
exponential decay of the kernel of S gets used. 

As already mentioned, except for the replacement of LemmaEJby Lemma|31 the proof of 
Proposition [21 consists of simply mimicking the discussion of the proof of Theorem ^ given 
in Section [21 



STEP 3. r Contains Projections onto GP Minimizers. 

Let ^'^P C L'^{M?) denote the set of all minimizers of the GP functional ©. We now 
consider the special case where S = — 1</>)(</>| for some G <^^^ . In this case, we claim that 



1 



hm j^El^'\N) = E^"- {a) - X 



(114) 



for any A > 0. Given Theorem^ the lower bound is trivial in this case. The upper bound 
can be derived in the same way as the upper bound for Theorem ^ in j23j . The arguments 
there also apply to this case, and the expectation value of S in the trial state can easily 
be estimated using the methods in [SI 115). (In the non-rotating case, this was carried out 
in [21].) 

Taking the derivative of ()114() at A > 0, Griffiths' argument [Hl^l implies that the one- 
particle density matrix of a ground state of H^^^i converges to \4>){'P\ as — > 00 in this case. 
Hence, by a similar 'diagonal' argument as at the end of the proof of part (i) of Theorem [21 

{A S) 

we can find a sequence A^ with A^ ^ as A — > 00 such that the ground state of H,^ ^ 
represents an approximate ground state sequence for the A = problem, and its reduced 
one-particle density matrix converges to |i;^>)(i;^>| as A ^ 00. This shows that \(f)){(f)\ G F for 
any (f) S 
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(Remark: The claim of this subsection can in principle be proved by simply constructing 
an appropriate approximate ground state. However, although the one-particle density matrix 
of the trial state used in [23, converges to \4>){4>\ as N ^ oo, this does not immediately imply 
that G T since the trial state is not symmetric! This explains the somewhat different 

reasoning in this subsection.) 

We note that also G Text for all (f> G <I>*~^^. This follows from the fact that all 

elements of T are positive operators, and a rank one operator cannot be written as a non- 
trivial convex combination of two positive operators. In the next subsection, we will show 
that all elements of Fp^t are of the form with 6 G <^^^. 



STEP 4. Proof of Parts (ii) and (iii). 

For a given 7 G F, let 7^^ be an approximate ground state sequence for H^^, with 7^^^ 7 
as ^ 00. By Proposition [21 we have that, for any bounded hermitian operator 5* and any 
A G M, 

i?«P(a) + ATr57 = Jim^^TV4^^)7^ >ii;gP)(a). (115) 
Upon dividing by A and letting A ^ 0, this yields 

Tr 57 > lim ^^^^ , — . (116) 

~ A\o A 

We claim that 

lim ^ = min . (117) 

Using (f> G ^"^^ as a trial function, we immediately see that E^^^^{a) < i?^^(a) -|- X{(f)\S\(j)) 
for all (f) G For the other direction, we use a minimizer of ^(^5) as a trial state for S'^^ . 

As A ^ 0, this sequence of minimizers will have a subsequence that converges strongly to 
a minimizer of <f^^. Hence, for some (/> G <I>*^^, limA\o A~-'^(i?^^'g,j (o) — i?^^(a)) > {(j)\S\cl)), 
which proves our claim. (Note that this argument also proves that the right side of H117|) is 
a true minimum and not merely an infimum.) 

We have thus shown that, for every bounded hermitian operator S, and every 7 G F, 

TV57 > min (0|5|(/>) . (118) 

Replacing S by —5*, this also implies that TrS'7 < max^g<^GP (0|5'|(/>). Inequality (|118|1 is the 
key to the proof of statements (ii) and (iii) in Theorem [21 

Let Pn be a rank n projection, and let P„F = {Pn^Pn ■ 7 G F}. When 7 is a bounded 
operator on L^(M'^), P^^Pn can be identified with an n x n complex matrix, and hence with 
a vector in M'^"^. We make this identification (denoted by l in the following) in order to 
be able to use finite-dimensional convexity theory (see, e.g., [Sn]). Note that l is linear and 
continuous, and hence the set i?„ = tPnT = {tPn^Pn : 7 G F} is a closed convex subset of 
M^" . An exposed point 20 of a convex set C C M™" is an extreme point p of C with the 
additional property that there is a tangent plane to C containing p but containing no other 
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point of C. (For an example of points that are extreme but not exposed, let C C be a 
square with each corner rounded off into a quarter of a circle. The extreme points are all 
the points on the four quarter-circles, including their endpoints, but the endpoints are not 
exposed.) 

An equivalent way to say this is that an exposed point p m C C is characterized by 
the existence of a vector a G (a normal to the tangent plane) such that 

{a,p)<{a,b) for ah 6 gC, (119) 

with equality if and only if b = p. (Here, (• , •) denotes the standard inner product in M*".) 

For a fixed n, an exposed point of Bn C M^"^ corresponds to some PnjPn £ Pn^- This 
density matrix 7 may not be unique and it may depend on n, but this is of no concern to 
us. 

We note that our space of density matrices is a complex space and, therefore, we have 
to translate H119|) to this setting. For any two bounded operators 7, 7' (not necessarily in 
r), the real inner product (• , •) becomes 

{iPn^Pn, iPniPn) = 3^ Tr(P„7tp„7') , (120) 

where 7^ is the adjoint of 7. Translated to our original space, this means that if Pn^Pn is 
an exposed point of PnF, then there exists an operator S (with PnSPn = S) such that 

3f? Tr S7 < 3ft Tr S7 for all 7 e F , (121) 

or, equivalently, there exists a hermitian S such that 

Tr S'7 < Tr ^7 for all 7 G F . (122) 

Note that, by definition, equality holds in H122j) if and only if P„7P„ = Pn'jPn- We now 
use inequality H122p . with 7 = where (f) € ^"^^ minimizes {(j)\S\(f>) among all GP 

minimizers. We know from Step 3 that this 7 is an element of F. The inequalities ()118() 
(applied to 7) and H122|) for this special choice of 7 together imply that there is actually 
equality in this case, and thus that PnjPn = Pn\(p){(p\Pn- That is, all exposed points of P^F 
are of the form Pn\4>){4>\Pn, with G . 

We can go further and conclude that all extreme points in P„F are of this form, not 
only the exposed points. This follows from the fact that the set of GP minimizers is closed, 
together with Straszewicz's Theorem [2{)\ Thm. 18.6] which states that the exposed points 
are a dense subset of the extreme points. 

Caratheodory's Theorem j2Ul Thm. 17.1] implies that every P„7P„ G PnF can be written 
as a convex combination of 2n^ + l extreme points. That is, there exist Aj > with Aj = 1 
such that 

(2n2+l \ 
M<l^i){(l)i\\ Pn, (123) 

with (pi G <I)^^ for all i. This equation defines an atomic (i.e., point) measure duni^P) 
supported on the (compact) space of projections onto GP minimizers. Let us provisionally 
call this space A, with the intention of showing that Fgxt = A. 
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For every ip with PnV' = "0 we have thus shown that 

iijhii^) = [ dfir^immi^ with [ df,n{(t>) = i. (124) 

J A J A 

To complete the proof of Theorem [51 we wish to take the hmit n — > cxd in (|124jl . We choose 
Pn in such a way that Pn converges strongly to the identity as n ^ oo. The sequence d^n 
has a subsequence that converges weakly to some measure dfi with J^dfi = 1 (see jH vol. 1, 
Thm. 12.7 and 12.10]). This implies that, for ip in a dense subset of L^(]R'^) (namely, those 
ip for which PnV' = foi' some n), 

(iPhli^) = [ dfi{4>)\{m\^ with [ dfi{(t>) = l. (125) 
J A J A 

Since H125() holds for a dense set of V', it actually holds for all ijj by continuity. That is, 
7 = J^dfi{(j))\(j)){(j)\ in the weak sense. 

Note that there is a representation 1)125(1 for 7 G Text (since there is such a representation 
for all 7 G r). It is not hard to see that for an extreme 7 the corresponding Borel measure 
dfj, must be an atomic measure at a single point in A. Another way to say this is that 
Text C A, which is exactly part (ii) of Theorem |21 (since we have already proved in Step 3 
that A C Text). 

Part (iii) of Theorem [21 follows from H125|) . together with part (ii). This completes the 
proof of Theorem [21 

We conclude with the direct proof of l|ll|) . which was promised just after the statement 
of Theorem 2. We start with (|123() and choose Pn to be the projection onto the largest 
n eigenvalues of 7, with n large enoug h so that Tr|7 - P„7P„| < We now denote 

Pn = P, 1 — Pn = Q and B = Ai|0i)(0i|. From ()123|) (and a little algebra) we learn that 
-f - B = - B)Q - QBP - PBQ. Thus Tr I7 - 5| < TV {\Q-fQ\ + \QBQ\ + 2\QBP\). 
Obviously, Tr Q7Q < e^/8 and since TrS = Tr7 = 1, we also have Tr \QBQ\ = Tr QBQ = 
Tr (1 — P)B = Tr (1 — P)j = Tr Q"fQ < The remaining term can be bounded, using 

Schwarz's inequality, by {Tr QBQy/^{TT PPP)^/^ < e/Vs. This proves (fTTIl . 

Appendix: Heat Kernel Estimates 

In this appendix we derive an upper bound on the heat kernel for a general Schrodinger 
operator. This bound will show, in particular, that for any s > and a large enough 
(depending on s) 

Tr|2;|"e"(^-^) < 00 (126) 

if V{x) > Ciln(|x|) — C2 for some constants Ci > and C2. This property was used in 
the proof of Theorem ^ (Actually, in the proof of Theorem ^ we used only the cases s = 2 
and s = 4 (see Step 4 of Sect. [^J because we assumed A = A x, but (|126() permits the 
inclusion of a magnetic field with polynomial growth of ^.) 
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Our bound on the heat kernel fohows an idea of Symanzik [26| . Using the Feynman-Kac 
formula for the integral kernel, we can write 



,y) = j exp ^ dsV{ojis))^ , (127) 



where d^x,y denotes the conditional Wiener measure for paths uj going from x to y in time 
a. By Jensen's inequality we have, for any given path uj, 

exp(- [ dsV{Lo{s))] < - [ ds exp{-aV{uj{s))) . (128) 



^ a 

Therefore (using Fubini's Theorem) 



a Jo ^ 



gSAg-ayg(a-s)A| (^^y)^ (129) 



To evaluate the trace in H126|) . we only need the heat kernel on the diagonal, i.e., for 
X = y. The integral kernel of e*^ is given by 

J.(x-.)-^^eH-^l^/-, (130) 

which leads to 

{e^^-V-)A}(,,,) = ^_l^ |^^d,e-^(^)exp(-|x-,|V4t) , (131) 

where t is defined byl/f = l/s + l/(a — s). 

Let us change the integration variable from s to and introduce the function 

2 /""/^ 1 
K{x) = - dt jt{x) . (132) 

a Jo Jl- At a 



Then the bound ((T^ yields 



e"(^-^)(x, x) < — — ^ (e-<^^ * h^) (x) , (133) 

with * denoting convolution. Note that J ha{x)dx = 1. It is easy to see that ha{x) ~ 
exp(— |xp/a) for large Hence, if V{x) increases logarithmically with we see that the 
diagonal of the heat kernel decays at least as " for large Thus, we can choose 

a large enough to ensure that H126|) is finite. 

For the proof of Theorem [21 it is necessary to extend this result to the case where — A + 
is replaced by —IS. + V+K, with K a finite rank operator. As explained there, we can restrict 
ourselves to the case when K has exponentially decaying eigenfunctions. I.e., we can assume 
that the kernel of which we denote by K{x,y), satisfies a bound 

i^(x,y) < Se-^d^l+ls'l) (134) 
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for some constants B > and D > 0. Again we want to show that, for any s > and a 
large enough (depending on s), 

Tr|x|^e"(^-^-^) < oo (135) 

if V{x) > Ci ln(|3;|) — C2 for some constants Ci > and C2. 

With the notation Lt = e*^^~^\ we can use the Dyson expansion to write 

^aiA-V-K) ^^^^ 5^(-l)" / dtodh ■■■dtn Lt,KLt,K ■ ■ ■ KLt^ . (136) 

n>l Jy:.U=a 



We have ah'eady derived an upper bound on the kernel of above. The kernel of the terms 
for n > 1 in the sum can be bounded as follows. First of all, the Feynman-Kac formula 
tells us that since F > we have the inequality Lt{x,y) < jt{x — y) for the kernel of Lj. 
Moreover, using (|134() and denoting by $ the function <l>(x) = V^e^^l^'l, we have 



{Lt,KLt,K---KLtJ{x,y)\<jt,*^{x) Wmum jt„ * ^{y) . (137) 



i=l 



Since Lt < I, we have |<I)) < ||^'||2- Denoting 

Ux) = Hulls' sup jt*^{x), (138) 

0<t<o 

we thus have 

\{Lt,KLt,K---KLtJ{x,y)\<U^)Uy)\\nf'- (139) 
The integral over the simplex in (|136() yields a factor a'^/nl, and hence 

e"(^-^-^)(x,y)| < e"(^-^)(x,y) + (e^H^H^ - l) U^)Uy) ■ (140) 

Since decays exponentially for large |x| this proves our claim (|135() . 
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